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Abstract: The purpose of this short review is to describe cosmological models with a linear 
inhomogeneous time-dependent equation of state (EoS) for the dark energy, when the dark 
fluid is coupled with dark matter. This may lead to a bounce cosmology. We consider 
equivalent descriptions in terms of the EoS parameters for an exponential, a power-law, or a 
double-exponential law for the scale factor a. Stability issues are discussed by considering 
small perturbations around the critical points for the bounce, in the early as well as in the 
late, universe. The latter part of the paper is concerned with dark energy coupled with dark 
matter in viscous fluid cosmology. We allow the bulk viscosity ( = ((H, t ) to be a function 
of the Hubble parameter and the time, and consider the Little Rip, the Pseudo Rip, and the 
bounce universe. Analytic expressions for characteristic properties of these cosmological 
models are obtained. 
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1. Introduction 

The current procedure for explaining the observed expansion for the universe [1,2] is to introduce 
a dark energy fluid with negative pressure and negative entropy [3-5]. According to the astronomical 
observations the dark energy currently accounts for about 73% of the total mass/energy of the universe 
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and only 27% of a combination of dark matter and baryonic matter [6]. The dark energy universe may 
have very interesting implications in the future [7]. A recent review of dark energy is given in Ref. [8]. 
A different way of accounting for the dark energy without any extra components is the modification of 
gravity (for recent reviews, see Ref. [9,10]). Note that in such a case there is the possibility of a unified 
description of early-time inflation and late-time acceleration as first realized in the Nojiri-Odintsov model 
[11,12] of F(R) gravity (for more models of unified inflation with dark energy in modified gravity, see 
Ref. [9]). In fact, the inhomogeneous fluid cosmology may also represent specific versions of F(R) 
gravity [11]. 

In the coupled phantom/fluid model, dark energy and dark matter are usually described by assuming 
an ideal (non-viscous) fluid with an unusual equation of state (EoS). New forms for the EoS for general 
(inhomogeneous/imperfect) dark fluid models were considered in Refs. [13-18]. The equation of state 
dark energy parameter w is known to be negative: 


Pd ^ n 
w — — <0, 

Pd 


( 1 ) 


where is the dark energy and pn is the dark pressure. According to present observational data the 
value of w is w = — l-CMi^io [19]. Various possible scenarios for the evolution of the universe discussed 
in the literature are concerned with the so-called Big Rip [20,21], the Little Rip [22-29], the Pseudo Rip 
[30] and the Quasi Rip [31] cosmologies. 

On the other hand in the early universe there is the possibility of a matter bounce cosmological model 
(see Refs. [32-34] and references therein). The universe goes from an accelerated collapse era to an 
expanding era over the bounce without displaying a singularity, what implies a cyclic universe. After the 
bounce the universe soon enters a phase of matter-dominated expansion. Recently the F(R) and F(T) 
gravity models in the presence of a bounce cosmology were discussed in Refs. [35-39]. We will in the 
following study similar bounce cosmologies. 

We will study a cosmological system with two coupled fluids: a dark energy component with a 
linear inhomogeneous EoS, and a dark matter component with a linear homogeneous EoS in a flat 
homogeneous and isotropic FRW universe. There exist other investigations dealing with the coupling 
between the inhomogeneous fluid dark energy and dark matter components [40,41]. We will explore 
bounce cosmological models where the scale factor is described by exponential as well as power-law 
forms and realize them in terms of the parameters of the EoS for the dark fluid. Then we will consider 
the stability of these models in the early universe by considering perturbations to the first order. Also 
we will consider a bounce cosmological model when the scale factor has a double exponential form. 
The bouncing behavior, both in the early universe and in the late-time universe, will be presented. The 
stability conditions for this model are discussed. 

For an introduction to bounce cosmology in general, the reader may also consult Ref. [42]. 

In Section 4 we study examples of Little Rip, Pseudo Rip, and bounce cosmology, described in 
flat Friedmann-Robertson-Walker spacetime when the cosmic fluid is viscous and is coupled with dark 
matter. We find corrections to the thermodynamical parameter, and to the bulk viscosity, in the equation 
of state for the dark energy. 
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2. Bounce cosmology via an inhomogeneous fluid 

In this section we construct bounce cosmological models induced by coupled dark energy in terms of 
the parameters in the equation of state. 

Let us consider a universe filled with two interacting ideal fluids: a dark energy component and a 
dark matter component in a spatially flat Friedmann-Robertson-Walker metric with scale factor a. The 
background equations are given by [7]: 

[ P + 3 H(p + p) = -Q, 

\ Pm + 3 H(p m + p m ) = Q, (2) 

[ H = - ! y(p + p + p m + Pm), 

where H = d/a is the Hubble rate and k 2 = 87 tG with G denoting Newton’s gravitational constant; p, p 
and p m , p rn are the pressure and the energy density of dark energy and dark matter correspondingly; 0 
is a function that accounts for the energy exchange between dark energy and dark matter. Here a dot 
denotes derivative with respect to cosmic time t. The Friedman equation for the Hubble rate is given by 
[7]: 

H 2 =j(p + p m ). (3) 

We will investigate the cases where the scale factor has an exponential, a power-law, or a double 
exponential form, and study bouncing behavior around t — 0 . 

2.1. Exponential model 

Let us first assume a bounce cosmological model where the scale factor a has an exponential form 
[35]: 

a = exp (at 2 ), (4) 

a being a positive constant. The instant t — 0 (a — 1) is taken to be the instant when bouncing occurs, 
in accordance with earlier treatments (cf., for instance, Ref. [35]). The Hubble parameter becomes 

H = 2 at. (5) 

We suppose that the dark energy component obeys an inhomogeneous EoS [13]: 

p = w{t)p + k{t). (6) 

The dark matter component is taken to obey a homogeneous EoS [ 8 ]: 

Pm = WPm■ ( 7 ) 

We will in the following assume that the dark matter is dust matter, p m = 0, corresponding to w = 0. 
Now consider the gravitational equation of motion for the dark matter: 


Pm "F 3 Hp m Q- 


(8) 
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We choose the interaction term between dark energy and dark matter in the following form, depending 
exponentially on time: 



(9) 


with Q o a constant. Thus Q(0) = Q 0 when t = 0. The solution of Eq. (8) becomes 



( 10 ) 


where p m { 0) = po■ According to this the dark matter energy increases linearly with t for low t, and fades 
away when t —y oo. 

Taking into account Eqs. (2, 3) we obtain the gravitational equation of motion 



( 11 ) 


As the time development of dark matter is known from Eq. (10), Eq. (11) can be regarded as the equation 
of motion for the dark energy. 

Let us choose the parameter w(t) in the form 


, , „ 6k 2 tt2 

w(t) = -1 - —H , 


( 12 ) 


with 5 a positive constant. Thus w(t) < 0 always, corresponding to the phantom region. When t —> 0, 
w(t) —> —1 (the cosmological constant case), while when t —> oo, w(t) —> — oo. 


Solving Eq. (11) with respect to A(t) (the time-dependent cosmological ’constant’) when w(t ) is 
taken to have the form (12), we obtain 



(13) 


Thus when t —> 0 one sees that A —>■ —4 a/k 2 — p 0 /2 (cf. Eq. (10)), while when t oo, A (t) —> 
5(2at) 4: —y oo. 

Concluding this subsection, we have explored the exponential model of bounce cosmology in terms 
of time-dependent parameters in the EoS, taking into account the interaction between dark energy and 
dark matter. 

2.2. Power-law model 

We next consider the case where the scale factor has the following form [35] (see also [41]): 



(14) 


where a ^ 0 is a constant, t is a reference time, and q = 2n, n 6 Z. 
In this case the Hubble parameter becomes 


a(g/t)(t/t) q 1 
a{t/t) q + 1 


(15) 
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In analogy to the case above, we assume also now that the dark matter pressure is zero, p m = 0. 
The derivative of H with respect to cosmic time is equal to 


t~\q- 1) _ H 
a(t/t) q + 1 q 


(16) 


Let us choose the interaction term between dark energy and dark matter in the form 

Q = _ Qe _ 

[a(t/t) q + l] 3 ' 


(17) 


As before, t = 0 means the bouncing time, corresponding to Q{t = 0) = Qq. When t —> oo, the 
interaction term Q(t ) —* 0. 

Solving the gravitational equation of motion (8) for the dark matter, we find 


Pm(f) 


P0 + Q 0 ^ 
[a(t/t)i + l] 3 ’ 


(18) 


where p m (0) = p 0 - 

Let us assume that the thermodynamic parameter w(t) in the EoS (6) has the same form (12) as 
before. Taking into account Eqs. (15)-(18) we obtain from Eq. (11) the following expression for the 
cosmological ’constant’: 


A (t) = 5H 4 


2 H r t~\q- 1) 

k 2 a(t/t) q + 1 



Pm- 


(19) 


Near t = 0 for n ^ 1 we have A —>• p 0 . In the case n = 1 we have A —* —4a/ k 2 P + p 0 . 

As a result we have explored how the power-law expression (14) for the scale factor, in addition to 
the form (17) for Q and the form (12) for w(t), influence the parameter A(t) in the EoS for the bouncing 
universe. 


2.3. Double exponential model 


Now we will consider a double exponential model where the scale factor has the form [35] 


a = exp(r) + exp(F 2 ), 


( 20 ) 


where Y = (t/t) 2 and t is as before a reference time. In this model there is a unification of the bouncing 
behavior both in the early universe and in the late-time cosmic acceleration. 

The Hubble parameter is equal to 


2 : exp(y) + 2Y exp(F 2 ) 

" t exp(Y') + exp(F 2 ) 


We suppose that the interaction term Q has the form 


Q 



( 22 ) 


with Q = Q 0 at t = 0. 
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The we can find the solution of the gravitational equation of motion for the dark matter [13]: 

Pmit) = ^ + QiVY, (23) 

a 6 

from which 

Pmit) = Q - 3Hp m (t). (24) 

If we take the parameter w(t) to have the same form (12) as before, we can solve Aft) from Eq. (11): 


A(t) = 5H* - - 


iW + (j 


Y 


1 + 


(1 + 4Y' 2 ) exp(Y' 2 ) 


+ wp n 


(25) 


When t —> 0 the cosmological ’constant’ A —>■ —2 /(k 2 t) — p 0 . 

As result we have constructed the double exponential model in the terms of the EoS parameters for 
the coupled dark energy. 


3. Stability of the solutions 


In this section we examine the stability of the solutions obtained in the previous section. First of 
all, we establish the formalism of this method. The evolution equations (2) can be expressed as a plane 
autonomous system [7] 


W = -[i+p'M] 
* = -[i -p'M] 


3 + 2 Hp k 

X 

3x + 2 Hp 

: y \ 


+ 3y[2x + (1 + w m )(1 - x - y)], 


= “ [ 3 * + 1(1 + ®m)(l - x ~ t/)] . 


here we have defined nondimensional quantities, 

k 2 p k k 2 Pv 

X ~ 3 H 2 ’ V ~ 3 H 2 ’ 

where p k and p v are the "kinetic" and "potential" terms, 

Pk=^{p + p), Pv = ^(p-p), 


(26) 


(27) 


(28) 


and N = in a, p'{p) = dp/dp. 

We shall investigate the stability of the critical points in the bounce models. Let us consider the early 
universe (1 —y 0), where bounce cosmology can be realized. We assume for simplicity that the bounce 
takes place at t — t*, so that t* < 0 (t* > 0) represents the contracting (expanding) phase. The bounce 
time /;* is assumed to be very small. In this model the universe has a collapsing era for t < A, and 
thereafter initiates an expanding era. 

With these assumptions the autonomous system (26) for the exponential, the power-law, and the 
double exponential models can be written on the common form 


% = 3x(l + x-y), 

*L = 3[ y (l+x-y)-2x\. 


(29) 
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Setting dx/dN = 0 and dy/dN = 0 in Eqs. (29), one obtains for the critical points: (x' c ,y'c) = (o,o) 
and (x", y") = (0,1). 

Consider now small perturbations u and v around the critical point ( x c .y c ), i.e., 


X = x c + u, y = y c + v. 


(30) 


With (x' c , y' c ) = (0, 0) the substitution of (30) into (29) leads in the linear approximation to the first order 
differential equations 


— = 3 ?/ 
dN 
dv 


% = 3(« - 2u). 

The general solution of (31) for the evolution of linear perturbations can be written as 


(31) 


u — Ci exp(31V), 
v = (c 2 + c 3 lV) exp(31V), 


(32) 


where ci, c 2 and c 3 are arbitrary constants. 

In the limit N —» oo, u, v —> oo. There occurs an asymptotic instability of the solution around this 
critical point. This instability may be a consequence of of the influence from dark matter on dark energy. 
In the case (x", y”) = (0,1) the autonomous system (29) in the linear approximation becomes 


The solution of (33) has the form 



0, 

—3 (u + v). 


u = Ci, 

v = c 2 + c 3 exp(-3!V). 


(33) 


(34) 


In the case c 3 = 0 we obtain a one-parametric family of points at rest on a straight line. The resting point 
(0,1) is stable. 

This concludes our study of stability of the solutions of the bounce cosmological models. 


4. Dark energy coupled with dark matter in viscous fluid cosmology 

Our second theme of the present review is to describe the theory of dark energy when it is coupled 
with dark matter in viscous cosmology. Cosmological models treating dark energy and dark matter 
as imperfect fluids with an unconventional EoS were considered in Refs. [13,14,18]. In these cases 
a viscous fluid is just one particular option. As is known, the case of a nonviscous (ideal) fluid is 
an idealized model, which is useful in many instances but in some situations incorrect. Viscosity, for 
example, is an essential ingredient when considering turbulence effects in cosmology [43]. From a 
hydrodynamic viewpoint the viscosity effect describes the deviation from thermodynamic equilibrium 
to the first order. The influence from bulk viscosity in the cosmic fluid plays an important role for the Big 
Rip singularity [13,44], or in type II, III and IV Rip singularities [13,45]. The Little Rip and Pseudo Rip 
phenomena and the bounce cosmology interpreted as an inhomogeneous nonviscous dark fluid coupled 
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with dark matter, are considered in Refs. [40,46]. Cosmological models in which the modification of 
gravity is described in terms of a viscous fluid are explored in Refs. [47] and [48]. 

A key reference for the material in the following is our recent paper [49]. 

4.1. Little Rip 

Among the theories of future singularities the most dramatic one is that describing the Big Bang 
scenario, in which the universe encounters a singularity in a finite time. There are also milder variants of 
the future singularity phenomenon, namely the so-called Little Rip and the Pseudo Rip. We now review 
these scenarios briefly. 

First consider the Little Rip. It is characterized by an energy density p increasing with time, but in 
such a way that an infinite time is required to reach the singularity. It corresponds to an equation-of-state 
parameter w < — 1, but w —* — 1 asymptotically. 

Let us take the Hubble parameter to increase exponentially with time, cf. [24], 

H = Hoe xt , Hq > 0, A > 0, (35) 

where H 0 = 77(0), t = 0 denoting present time. 

Assume the dark matter to be dust, p m = 0, so that the gravitational equation for dark matter reduces 
to 

pm + 3 Hp m = Q. (36) 

With the same coupling as in [9], 

Q = 5Hp m , (37) 

the solution of Eq. (36) for dark matter becomes 

-Hj , (38) 

where p 0 is an integration constant. If 5 is small (< 3), p m —> 0 when t —> oo. 

One can choose different forms for the equation of state for the viscous fluid. We adopt the 
inhomogeneous form; cf. [47] 

p = w(p)p-3HaH), (39) 

implying that the bulk viscosity depends on the Hubble parameter. Further, we choose w(p) to have the 
form; cf. [47] 

w(p) = A 0 p“ _1 - 1, (40) 

where A 0 0 and a > 1 are constants. This leads to the the gravitational equation of motion for dark 
energy 

6 77 3 

-p- - Pm + 3 H[A 0 p a - 3H((H)} = —6Hp m . (41) 

From this we obtain an explicit expression for the Hubble-dependent bulk viscosity 

C( H ) = ^ j^(Pm + A 0 p a ), 


pm{t ) = Poexp 


5- 


(42) 




Entropy 2015 , xx 


9 


where the second term contains the contribution from the coupling. 

When t —> oo, H —> oo exponentially, as does the density p. The second term to the right in Eq. (42) 
goes to infinity in this limit. 

It is of interest to consider also another variant of the Little Rip model where the Hubble parameter 
increases with time according to a double exponential, 

H = H 0 exp(Ce xt ), (43) 


where H 0 , C and A are positive constants. 

We take the interaction term Q to have the same form (37) as before. Solving the gravitational 
equation of motion (36) for dark matter, we find 


Pm{t') 


P o exp 



5-3 


A 


-H 


(44) 


Here p 0 is an integration constant. 

We keep the same form (40) for the dark energy thermodynamical parameter w(p) in Eq. (39) as 
before. Then we obtain from the gravitational equation of motion (41) for dark energy 


C (H,t) 


‘ZCX w A 0 p a p m 
3 k 2 3H 3H 


5 _ 

3 


(1 - e c ) + e c+xt 


(45) 


In this case the cosmic time occurs explicitly in the viscosity, not only indirectly via p, p rn and H. 


4.2. Pseudo Rip 

This is another soft variant of the future singularity scenario in which the Hubble parameter tends to 
a cosmological constant as t —» oo. That means, the universe approaches a de Sitter space. We will now 
analyze this model in analogy to the analysis of the Little Rip model. 


Assume first that the Hubble parameter has the form; cf. [24] 

H = H 0 -H 1 exp (—At), 


(46) 


where H 0 , Hi and A are positive constants, // (l > II\.t> 0. Thus in the late-time universe H —> 7/ 0 . 

With the interaction between dark energy and dark matter in the form (37), we find as solution of the 
gravitational equation (36) 


Pm(t) = Po exp 


(5-3) 





(47) 


with po an integration constant. 

If the thermodynamical parameter w(p) has the form (8) we obtain from Eq. (41) the bulk viscosity 


as 




2A 

W 2 



+ 


1 

3 H 


(Pm + Ao p“), 


(48) 


showing again the coupling. 
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Let us also consider, as a second example, a cosmological model where the Hubble parameter has the 
form; cf. [24] 


H = 


*L 
V3 


1 



VSAt 
2 x f 


(49) 


Here x 0 = sfpo corresponds to the present energy density, Xf is a finite, and A is a positive constant. 
In the late universe H —y Xf/y/ 3, and the expression (49) tends asymptotically to the de Sitter solution. 
When t —y 0, we obtain H —y x 0 /-s/3- 

For the energy density of dark matter we find 


Pm{t) = Poexp 




t + 


2 Xf 

y/ZA 



with p 0 an integration constant. 

In this case we derive the following bulk viscosity: 




A 

3 k 2 H 



+ 


1 

3 H 


(. Pm + A 0 p a ). 


(50) 


(51) 


The interaction between dark energy and dark matter implies corrections in the bulk viscosity. This 
contrasts the case where only dark energy was involved. 


5. Conclusion 

In the present paper we investigated bounce cosmological models in which we took into account the 
interaction between dark energy and dark matter. In terms of parameters of the equation of state, w(t) 
and A (t) (cf. Eq. (6)), we described bounce cosmologies when the scale factor is expressed by either 
an exponential, a power-law, or a double exponential form. In all cases, we adopted the form (12) for 
the parameter w(t). It turned out that near the bouncing instant t — 0, the expressions for A (t) reduced 
essentially to the initial density of dark matter, p 0 . 

Not very much is known from observations about the form of the interaction term Q. The analytic 
forms given above, in Eqs. (9), (17), and (22), were motivated chiefly by mathematical tractability. 
Physically, we have as checking points that Q —> 0 when t oc. 

We also analyzed in a linear approximation the stability of the stationary points against perturbations 
and showed the existence of a stable point, and attractor solutions, for these models. 

In the second part of our paper we studied examples of Little Rip, Pseudo Rip, and bounce cosmology, 
described in flat Friedmann-Robertson-Walker spacetime when the cosmic fluid is viscous and is coupled 
with dark matter. We found corrections in the thermodynamical parameter, and to the bulk viscosity, in 
the equation of state for the dark energy. 

It could be mentioned that there are common properties of the corrections present in these models. 
Thus A Cm(H, t ) = p m /(3 H) represents the correction to the bulk viscosity from the coupling. Further, 
A C P (H, t ) = ( A 0 /3H)p a is related to the choice of thermodynamical parameter w(p) in the equation of 


state. 
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From a general perspective the present theory may be considered as inhomogeneous cosmic fluid 
theory with changeable thermodynamical parameter w coupled to dark matter [50], applied to the Little 
Rip, the Pseudo Rip, and the bounce phenomena. 

As is known, cf. [9], a viscous fluid may be understood also as a modified gravity model, for instance 
of the F(R) type. It is moreover known that F(R) gravity may provide a unification of early-time 
inflation with a special version of dark energy, as was proposed in Ref. [12], or with the nonlinear model 
proposed in Ref. [51]. Having that in mind, we expect that there is a natural possibility to unify these 
epochs with inflation in an extended viscous model. 

There are different types of bounce universes. We have presented some examples that are in 
conformity with our model. Other types of bounce universes, such as those encountered in Matter 
Bounce Loop Quantum Cosmology for instance, will be discussed elsewhere. 
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